Critical indices of random planar electrical networks 
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We propose a new method to estimate the critical index t for strength of networks of random 
fused conductors. It relies on a recently introduced expression for their yield strength. We confirm 
the results using finite size scaling. To pursue this commitment, we systematically study different 
damage modalities of conducting networks inducing variations on the behavior of their nonlinear 
strength reduction. 
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I. INTRODUCTION 

Random resistor networks are used to model a vari- 
ety of physical phenomena rangingfrom the properties 
of inhomogeneous media P,[3,UI1,[5L to metal insulator 
transitions [|| , dielectric breakdown 3, HI, [t| [To| . the role 
of percolation in weak and strong disorder ll ll. Il2l . [l3| , 
and the strength of trabecular bone [3, EH, Tla ] . Essen- 
tially, this study intends to develop a method to establish 
a relationship between the mean strength of bone and its 
mass. In fact, large bones consist of an outer compact 
shaft and an inner porous region, i.e., a trabecular archi- 
tecture whose structure is reminiscent of a complex sys- 
tem of disordered networks. The bone strength depends 
on many factors, i.e., architectural characteristic (level of 
connectivity), perforation, thinning, anisotropy, as well 
as subarchitectural properties of bone like the mineral 
content, the density of the diffuse damage. A significant 
mechanism for loss of trabecular mass is through the re- 
moval of individual struts, due to traumatic events. Me- 
chanical studies on ex vivo bone samples have shown that 
trabecular networks from patients with a broad range of 
age fracture at a fixed level of strain [13], even though 
the corresponding fracture stresses exhibit large varia- 
tions. These observations have motivated the induced 
fracture criteria utilised in our models. 

We consider the lattice network of conducting disor- 
dered elements to study the essential non-linear and irre- 
versible properties of the electrical breakdown strength. 
Subjecting the system under extreme perturbation, its 
electrical properties tend to get destabilised and so that 
failure breakdown occurrences. In fact, these instabili- 
ties in the system often nucleates around disorder, which 
plays a major role in the breakdown properties of the 
system. The growth of these nucleating centers, in turn, 
depends on various statistical properties of the disorder, 
namely the scaling properties of percolation structures, 
its fractal dimensions, etc. By increasing the percentage 
number of the network nonconducting elements, its con- 
ductivity decreases, so also does the breaking strength of 
the material, the fuse current of the network decreases on 
the average with the increased concentration of random 
impurities. 



We are considering here the problem associated with 
failures in disordered system under the influence of an 
electrical field. In spite of this framework is shown sim- 
pler than the one related to mechanical failure of frac- 
tures, however all these cases of failures present some 
common features. In particular, certain features near 
breakdown agree with those of resistor networks close to 
the percolation threshold [H, El HJ EH- This article 
presents a systematic new method to estimate critical 
indices by quantifying the reduction of current in fused 
resistor networks. 

The model and the expression which relate the 
strength reduction to the statistical properties of disor- 
der network system is set up in section II. We focus on 
fractured conducting networks to examine many possi- 
bilities for. Initially, it is supposed that disorder just 
arises from random percolation. Both cases are consid- 
ered, the isotropic and the anisotropic conducting net- 
works. Next, it is also considered that the failure current 
of the conducting network are not the same, but it is be- 
ing uniformly distributed in a range, in addition to ran- 
dom uniform percolation. This issues are presented and 
quantified in Section III. Finally, Section IV is devoted 
to conclusions. 



II. THE MODEL 

We take square networks consisting of fused conduc- 
tors that fail when the potential difference across them 
reaches a pre-set value; i.e., the breakdown current of 
an element is proportional to its conductance. Typi- 
cally, failure of an element increases currents on neigh- 
boring conductors, enhancing the likelihood of their fail- 
ure [23, [H| . We study the yield point at which the exter- 
nal current initiates the first failure. The peak currents 
on a network show similar behavior [24j . 

Consider first, a complete square network of size M x 
M, with the top and the bottom edges at potentials Vo 
and respectively. Assume that each electrical element 
in the network fails when the potential difference across 
it reaches a value Vb- We can calculate the current 1(0) 
flowing through the network using Kirchhoff's laws. As 
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Vb increases, currents through the conductors, as well as 
7(0), will increase until the yield point 1(0) — l y i e id(0), 
where the first failure occurs. 

Next consider a network, where we have attempted to 
remove elements with a probability p. Denote the yield 
current of such a network by Iyieid (p) ■ Typically, I V i e id (p) 
decreases with increasing p, and vanishes as (po ~~ pY 
when p approaches the percolation threshold po (=0.5 
for isotropic square networks), t is the critical index, 
with reported values between 1.1 and 1.43 under differ- 
ent scenarios of dama ge a nd symmetries of the network 

[H Hi M, m m MM m \m ■ 

There have been several proposals related to the form 
of the reduction of strength of a network due to a ran- 
dom removal of element [3 EH- Here we test a recently 
proposed expression for I yie id{p) HI: 



r(p) = 7 » ieW ^ 



1 



I V ield(0) 1 + ai z t/2 + a 2 z l ' 



(1) 



where z = log (AT)/ log(^). Here, N(= M 2 ) is the num- 
ber of nodes in the original network, po is the bond per- 
colation threshold for the class of network considered, 
and a\ and a 2 are constants that depend on model pa- 
rameters as discussed below. Observe that, as p — » po, 
the yield strength, r(p) — » (j>o — p) 1 . We conjecture that 
Equation (1) is valid throughout the range p G [0,po), we 
use this conjecture to estimate po and the critical index 
t ; then validate it using finite size scaling method. 

We compute the yield current for a given network as 
follows: given the conductance of all electrical ele- 
ments in the network and the potential V of the top 
layer of nodes, we use Kirchhoff 's laws to determine the 
currents ik through each element and the potential differ- 
ences Vk across them. We denote the largest of the latter 
by V m ax- The current I(p) passing through the network 
is the average of all currents through electrical elements 
belonging to a fixed horizontal layer. The yield current 

IS Iyieldip) = I(P) X V b /V rnax . 



III. RESULTS AND DISCUSSIONS 

Our computations used networks with sizes M = 
20, 30, 40, 50, 70, 80, 90, 120, 140 and 160, and sample 
sizes of 10, 000 for M = 20, 30, 40, 50, sample size of 
2, 500 for M = 70, 80, 90 and finally sample size of 1, 000 
for M = 120, 140, 160. Figure Q] shows the behavior of 
r(p) for the 160 x 160 networks, where the error bars 
show standard errors. Although fluctuations St in r(p) 
decrease as p — > po, the relative fluctuations (St/t) in- 
crease. The solid line shown in Figure [1] represents the 
best fit to Equation (1) with parameters a%, Po and 
t in a 160 2 conducting network. To determine their val- 
ues, a x = -0.1043 ± 0.005, a 2 = 0.061 ± 0.003, p = 
0.513 ± 0.002 and t = 1.228 ± 0.015, we used the 
Laven berg -Marquadt method to implement the nonlin- 
ear fit [35j . We must fit t and po because they depend on 
the size of the network. 




FIG. 1: Strength reduction in a 160 2 network due to a random 
isotropic removal of a fraction p of conductances, averaged 
over one thousand trials. Numerical results are shown by 
boxes along with the statistical error, while those obtained 
by fitting to expression (1) are shown as a continuous line. 



Next, we determine how po and t change with the net- 
work size. We express the parameters as a function of 
x = M^ 1 / 1 ' where v (=4/3 for 2D square networks) is 
the universal correlation length exponent [3ol l37| ; i. e, x 
is the inverse of the mean size of the largest domain in 




FIG. 2: Variation of parameters t and po versus the inverse 
correlation length (x = Af -0 ' 75 ), for networks of size M = 
20 to M = 160. All elements in the initial networks have 
conductance of 1 unit. Note that, both critical values for 
t(= 1.066 ± 0.056) and p (= 0.500 ± 0.006) are getting in 
the limit when x — * 0. 
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FIG. 3: The left side of the figure shows the scaling function 
(r) behavior as a function of the network correlation length 
(£) , and its x error as a function of t for fixed unit conduc- 



tances in the network (right side), 
gives the critical exponent. 



The minimum of \ i s t, 



a network of size M x M. Figure shows the values of 
t(x) and po(x) , along with the error estimates. 

We estimate the value of i(0) corresponding to an in- 
finite network by first approximating t(x) by a rational 
function f(x)/g(x) (where f(x) and g(x) are polynomials 
of order 3 and 2 respectively) and extrapolating to x = 0. 
The values of the extrapolation corresponding to Figure^ 
are p (0) = 0.500 ± 0.006, t(0) = 1.066 ± 0.056, a x (0) = 
-0.106 ± 0.015 and a 2 (0) = 0.060 ± 0.004. The error 
estimate includes both errors at each M and those due 
to the extrapolation [35| . 

We now validate our results using finite size scaling 



to independently estimate i(0). According to the fi- 
nite size scaling ansatz, for the correct t(0), the rela- 
tionship between the rescaled variables r = M'(°)/ i/ x r 
and £ = M~ 1 / v {j>q — p) is independent of the system 
size M 38]. Although the finite size scaling ansatz need 
hold only for p — > pq (and < M~ x l v < 1), we find 
that the data collapses to a scaling function f (£) over 
the entire range p G [0,po]; see Figure [3] (left side). We 
then determine the best value for t(0): for any chosen 
value of i(0), we approximate the scaling function r(£) 
by a rational function f(x)/g(x) (where f(x) and g(x) 
are polynomials of order 3 and 2 respectively) and esti- 
mate the deviation of the data (Cfc,Tfc) from the scaling 
function by x 2 = ^2 k (jk — ^(Cft)) 2 - Here, the sum is 
over all available networks. Figure [3] (right side) shows 
the x 2 as a function of t(0); the best estimate, which 
we assume minimizes % 2 , is tpss = 1-07 ± 0.10, where 
the error estimate corresponds to doubling the % 2 value. 
This estimate agrees with that value we obtained from 
Equation (1). We tested all remaining estimates for pa- 
rameters using Equation (1) and the finite size scaling 
ansatz. 

Next we consider networks from which we removed el- 
ements anisotropically. We begin with a square network 
of unit conductance and remove elements in the horizon- 
tal and vertical directions with probabilities p^ = ftp 
and p v = p. Analysis of such networks for — 2.0 us- 
ing Equation (1) gives p = 0.3383 ± 0003, t = 1.3303 ± 
0.0189. Similarly, for = 1/2 , p = 0.6798 ± 0.0099, t = 
1.0396 ± 0.0553 (See Figure SJ). Finite size scaling for 
the two cases estimates give tpss = 1.3 ± 0.10 and 
tpss = 1-05 ± 0.10 respectively, in good agreement 
with Equation (1). The estimates for po f° r multiple 
/3's, shown in Figure [H agree with theoretical results for 
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FIG. 4: Values of t(/3) and po{/3) with f3, for anisotropic re- 
moval of conductance. Elements in the vertical and horizontal 
directions are removed with probabilities p and /3p. 




FIG. 5: Behavior of t(s) andpo(e) with e, when conductances 
of the initial network are chosen randomly within (1 — e, 1 + e). 
Bonds removal is isotropic (po remains at 0.5). 
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anisotropic networks [lj. The value of the critical ex- 
ponent changes little over the range from > 1.0 to 
P < 2.0. 

Now, we consider the slow breaking network process 
which consist of two subsequent independent random 
processes, i.e., after a random remotion of conductors; it 
is imposed to the network remaining elements a conduc- 
tivity value choosing at random in a range (1 — e, 1 + e); 
where e € (0,1). Note that, the isotropic case is re- 
trieved for e = 0. We find that the critical index de- 
pends on the value of e. We conduct our analysis for 
e = 0.00,0.25, 0.50, 0.75, 0.85, 0.95 and 1.0, and at each 
e, for network sizes considered earlier. Although the 
value of the critical point po is independent of e (Fig- 
ure [5]) , and the critical exponent increases with e (for 
e > 0.5) (Figure [5]), consistent with results for finite 
size scaling. 

Next, we consider networks whose remaining electri- 
cal elements randomly degrade as we remove conduc- 
tances. This problem relate to degradation of porous 
bone with aging [39| . Specifically, we consider networks 
whose conductance and breakdown decrease by a factor 
(1 — ap), where as before, p is the probability for an el- 
ement to be removed from the network. Once again, as 
earlier, we analyzed network sizes ranging from M = 20 
to M = 160 (Figure El). For a = 1.0 we find that p = 
0.509 ± 0.005, t = 1.287 ± 0.034, a% = -0.123 ± 0.006 
and a 2 = 0.067 ± 0.003. The estimated critical exponent 
using finite size scaling is tpss = 1-29 ± 0.10. 

Let us now stress a point on the essential of nonlinear 
and irreversible properties of the breakdown process. In 
Figure [7] are shown the strength reduction behavior orig- 
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FIG. 6: Variations in parameters t(a) and Po(a). When ele- 
ments are removed from the network with probability p, the 
conductance of those remaining are reduced simultaneously 
by a factor (1 — ap). The critical fraction removed elements 
remain 0.5. 
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FIG. 7: The deviation from linearity of the strength reduc- 
tion due to the removal of a fraction of conductance in a 
160 2 network. The right side shows the anisotropic removal 
of conductance with probabilities p and 3.0 p in the vertical 
and horizontal direction. The left side corresponds the case 
where the elements are randomly removed with probability p 
meanwhile the ones remaining are simultaneously reduced by 
a factor (1 - 0.75 p). 



inated by two different types of damage modalities. As 
can be seen, both deviates from linearity above the per- 
colation threshold in agreement with [4(| . The right side 
of the figure shows the strength reduction behavior for 
the anisotropic case for f3 — 3.0. Left side of the figure 
considers an ensemble of conductor which are randomly 
removed from the network while their remaining elements 
are diminished at random by a factor (1 — 0.75 p). For 
the slow breaking process, the strength reduction shows 
a similar behavior. 



IV. CONCLUSIONS 

We have computed critical indices for several classes 
of square networks of conductances using finite size scal- 
ing and Equation (1). For the isotropic case, t = 
1.066 ± 0.056 in agreement with the value t = 1.1 re- 
ported by Kirkpatrick (27|, Straley HH, HI], and Stinch- 
combe and Watson However, it is different from the 
values (close to 1.3) obtained from real space renormal- 
ization group methods [3l|, [H, EH • (This discrepancy has 
already been discussed by Straley [HI.) For anisotropic 
networks, we compare to an analysis of experimental and 
computational results by Han, Lee and Lee [46j and by 
Smith and Lobb |47|- Both groups found that t = 1.3 
when p = 0.33. Figure [H shows that for p = 0.33, the 
parameter j3 is 1.94 and t = 1.31. Further, the values of 
Po(/3) for anisotropic removal of conductances (Figure |4j 
is consistent with theoretical results of Redner and Stan- 
ley Q. 

Here, we have presented substantial evidence that criti- 
cal indices depend on the type of initial network and dam- 
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age modalities of conductances and the network. Indeed 
the critical exponent might be a convoluted exponent be- 
cause of two independent random processes are affecting 
the strength reduction of networks. Results from previ- 
ous studies have been shown to be isolated examples of 
our more general analysis of this problem. It would be of 
interest to develop a renormalization group based analy- 
sis to describe these more general damage processes. We 
have, thus far, not been successful in identifying such a 
theory. 
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